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Non-magnetic impurities are known as strong pair breakers in superconductors with pure d-wave 
pairing symmetry. Here we discuss d-wave states under the combined influence of impurities and 
competing instabilities, such as pairing in a secondary channel as well as lattice symmetry break- 
ing. Using the self-consistent T-matrix formalism, we show that disorder can strongly modify the 
competition between different pairing states. For a d-wave superconductor in the presence of a 
subdominant local attraction, Anderson's theorem implies that disorder always generates an s-wave 
component in the gap at sufficiently low temperature, even if a pure d^2_y2 order parameter charac- 
terizes the clean system. In contrast, disorder is always detrimental to an additional d^y component. 
This qualitative difference suggests that disorder can be used to discriminate among different mixed- 
gap structures in high-temperature superconductors. We also investigate superconducting phases 
with lattice symmetry breaking in the form of bond order, and show that the addition of impurities 
quickly leads to the restoration of translation invariance. Our results highlight the importance of 
controlling disorder for the observation of competing order parameters in cuprates. 



I. INTRODUCTION 

Although a consensus on the basic d^2_j^2 symmetry 
of the superconducting order parameter in most cuprate 
materials is now established 0], the possibility of fur- 
ther symmetry breaking due to the presence of competing 
ground states has been widely discussed, both on theo- 
retical and experimental grounds. 

One of the possible additional order parameters is a 
secondary superconducting gap of different symmetry. 
In the resulting mixed-gap state, the minimum free en- 
ergy is usually taken for complex order parameters like 
s-\-idx2_y2 and dxy + id^^-y^^ which implies the violation 
of time-reversal invariance (that we will call T-breaking), 
in addition to the broken U(l) symmetry of the super- 
conducting state. Experimentally, the observed splitting 
of the zero-bias conductance peak of tunneling exper- 
iments in the overdoped regime of YBa2Cu307_5 Q has 
been interpreted as evidence for a T-breaking mixed-gap 
state. However, the precise nature of this mixed gap is 
not entirely clear, as a number of studies reduce to ten- 
tative fits of the tunneling spectra using the different gap 
structures allowed by symmetry |5i . We also point out 
that the signature of such a mixed gap has remained elu- 
sive in other physical probes. Surface effects may play 
a role, as a small secondary gap can be induced by the 
proximity to interfaces 0, 13 j although other mechanisms 
based on competing pairing attractions |^ , magnetic 
fields effects ^3 and presence of magnetic impurities [TTI | 
have been proposed. It is worthwile mentioning that a 
quantum phase transition, associated with the passage 
from the pure (i^2_j,2-wave phase to the mixed supercon- 
ducting state, would provide a source of anomalous scat- 
tering of the nodal quasiparticles in the quantum critical 
regime IT^ . Photoemission measurements in opti- 
mally doped Bi2Sr2CaCu208+5 131 have indicated the 
presence of such critical scattering, however, experimen- 
tal details have not been fully sorted out. 

In electron-doped cuprates the situation appears to be 



different, as several experiments on Pr2-xCea;Cu04 indi- 
cated a rather sharp transition from a (i^ 2_., 2-wave to a 
s-wave state as a function of doping [ij, [ig. Therefore 
a mixed-gap phase could at best be present in a narrow 
region of the phase diagram. 

Of course, in addition other symmetries could be bro- 
ken in those materials, such as the lattice symmetry 
(denoted C-breaking in the following). Indeed, theoret- 
ical studies of models for doped Mott insulators have 
demonstrated the possibility of several different order- 
ing patterns, such as bond (or spin- Peierls) order j^, 
stripes 0, 0| , and staggered flux states 0| . On the 
experimental side, neutron scattering measurements re- 
vealed stripe formation in La2-xSra;Cu04 |23|. More 
recently, signatures of charge order were observed by 
tunneling experiments in the superconducting phase of 
Bi2Sr2CaCu208+5 jllEl IH, giving further evidence 
that competing orders may be generic to these materials. 

In this paper, we argue that disorder plays a significant 
role for competing instabilities in unconventional super- 
conductors, by possibly suppressing or enhancing the ten- 
dency for the system to form complex symmetry breaking 
patterns. A basic remark is that non-magnetic impuri- 
ties are known to strongly suppress the pure d^2_y2-wave 
state, while usual s-wave superconductors are extremely 
robust to the addition of static disorder, according to 
the famous Anderson's theorem 24]. This has imme- 
diate consequences for the competition of pairing states 
with different symmetry. First, a putative phase tran- 
sition such as d^2_j,2 — > s -|- idx2_y2 could be tuned by 
varying the disorder strength (e.g. using irradiation or 
addition of impurities). Second, we can expect that order 
parameters of s + idx2 _y2 and d^y + idx2_y2 types should 
behave quite differently to the addition of non-magnetic 
impurities, and this effect (which is discussed in more de- 
tail below) could be used to experimentally discriminate 
among the several gap structures discussed for cuprate 
and other unconventional superconductors. 

To gain some physical understanding, the first part of 
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this paper is devoted to the simplest theory of disorder in 
superconductors with a two-component gap. We argue 
that the BCS formahsm together with a self-consistent 
T-matrix approximation 25] is enough to capture the 
essential features of the problem, although quantum crit- 
ical dynamics 12], localization phenomena |2(tI | and non- 
uniform pairing patterns j27l | are not described by our 
calculation. In particular, some of our results obtained 
in the strong disorder regime should be taken with some 
caution due to the development of strongly inhomoge- 
neous regions [23|. Nevertheless, we will be able to de- 
termine the phase diagram as a function of the impurity 
concentration and of the ratio of two competing pairing 
attractions, which can be expected to be qualitatively 
correct. Previous studies of this problem 291 were de- 
voted to the calculation of the critical temperature in the 
Ginzburg-Landau formalism, and to our knowledge the 
solution of the BCS problem with disorder for s-|-«c?^2_j,2 
and dxy+idj;2_y2 superconductors was not investigated in 
detail. We shall demonstrate that disorder generates and 
stabilizes an s-wave part in the s-\-idx2_y2 situation. This 
result is due to the combination of two remarkable prop- 
erties of superconductors: the insensitivity of the s-wave 
gap to disorder (Anderson's theorem) and the presence 
of a finite density of states in a disordered dx2-y2-wsLYe 
state, leading inevitably to a secondary BCS instabil- 
ity in the s-channel at low temperature (a similar result 
was discussed in [s^, but here we are able to provide a 
more transparent physical picture). Because Anderson's 
theorem does not hold for a dxy + idx2_y2 superconduc- 
tor, the effect of disorder in this case is quite similar to 
the pure fia;2_j,2-wave case: superconductivity is rapidly 
destroyed, and no emergent order parameter can follow 
from impurity addition. Because of these two different 
situations, we propose that disorder can be used in prac- 
tice to distinguish the two debated scenarios in overdoped 
YBa2Cu307_5, and hope to stimulate new experimental 
studies. 

The second part of the paper concerns the sensitivity 
to disorder of superconductors with spontaneous lattice 
symmetry breaking, motivated both by stripes physics 
and possible charge ordering in cuprates. For simplicity, 
and because of the local nature of the T-matrix approxi- 
mation, we restrict our attention to superconductors with 
additional bond order, and demonstrate that the non- 
magnetic impurities have a general tendency to restore 
in an extremely fast manner the global translation sym- 
metry. 

The remainder of the paper is organized as follows: 
In Sec. ^ we concentrate on translation-invariant dirty 
superconductors with mixed gap symmetry, and de- 
velop the basic formalism, referring the reader to Ap- 
pendix ^ for technical details. We investigate in turn 
the s + id,j.2_y2 and d^y + id.j;2__y2 cases, discussing the 
evolution of the superconducting gap as a function of 
the impurity concentration and of the ratio of attrac- 
tive interactions in the two channels. In Sec. IIIII we 
study the more general problem of the influence of dis- 



order on superconductors that combine U(l) and lattice 
symmetry breaking. The complete derivation of the T- 
matrix formalism is given in Appendix ^ A general dis- 
cussion of disorder effects in cuprates is given in Sec. 
IV, together with a summary which will close the pa- 
per. Most of the numerical calculations are performed 
assuming strong impurities, i.e., unitary-limit scatterers. 
However, we expect the qualitative results to be insensi- 
tive to details of the disorder potential. 



II. T-BREAKING IN DISORDERED 
SUPERCONDUCTORS 

A. General formalism 

In this Section we will investigate the instabilities of a 
two-component superconductor using the weak-coupling 
BCS theory, owing to the fact that strong-coupling ef- 
fects would only modify quantitatively the results, ac- 
cording to previous results in the Eliashberg framework 
|3lj| . However, the quantum-critical dynamics near the 
phase transition between states of different pairing sym- 
metry |l2lll3| cannot be accounted for in our formalism, 
as fluctuations of the order parameter are discarded in 
the mean-field approximation. 

Introducing a general pairing attaction Vkk' < in the 
singlet channel, which corresponds to the Hamiltonian 

= E 'k^L^k. + E Vkk'i'l^i'Ui^j.k^i^k'^ : (1) 

k(7 kk' 

the complex superconducting gap A'^, + iA'^ obeys the 
BCS equation, 

A', + ^A'l ^~J2^kk'^J2 [{r' + iT^m\ ^c.„)] , 

k' n 

(2) 

where f3 is the inverse temperature and ujn — (2n-|- l)7r//3 
the Matsubara frequency. In the following, we will as- 
sign the two real gap functions AJ., A'^' to two different 
representations of the lattice symmetry group, keeping 
the global phase of the order parameter fixed. Further, 
g{k, iujn) is the single-particle Green's function of the 
Nambu spinor 'i'k = (V'feT' V'lfe|)- It can be decomposed 
as g = 5iT', with t', i = 1,2,3 being the Pauli matri- 
ces, and r'^ is the identity matrix. The Green's function 
is determined by Dyson's equation: 



i>2 + A;2^Af +?2 

where cj- is the quasiparticle dispersion and the contribu- 
tion of impurity scattering to the electronic self-energy is 
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S(fc,zci;„) = iw„)T'. We have precedently intro- 

duced the frequency-dependent renormahzed quantities: 



A'j, ^ Ei(iw„) , 
A','-t-S2(zc^„), 

Cfc -I- T^zilL^n) ■ 



(4) 
(5) 
(6) 
(7) 



We will determine S in the self-consistent T-matrix 
approximation, which employs the exact solution of the 
single impurity problem. If we note by v^w the static 
potential generated by a given impurity, the T-matrix 
obeys the Lippmann-Schwinger equation 25] : 

r(fc, fc', iw„) = Wfcfc'T^-f^ Vkk"T^g{k", iuJn)T{k", k' , «w„). 
k" 

(8) 

In the case of a finite impurity concentration, rii^p, the 
T-matrix approximation amounts to the neglect of inter- 
ference effects between different impurities. The main 
limitation is that localization effects [s^ are ignored, 
which may be important in low-dimensional systems (see 
[2^ and references therein); also, the occurence of inho- 
mogeneous states due to disorder cannot be accounted 
for [231 . Within the T-matrix approach, disorder aver- 
aging effectively restores translation invariance, and the 
disorder-induced self-energy is given by the relation: 



where (j) is the angle specifying the direction of k in polar 
coordinates. This results in s + id^2_y2 superconducting 
order: 



A. = A,, 



iAd cos(20). 



(13) 



Ordered states with a relative phase between the two 
components different from 7r/2, such as s -|- dx^-y^, are 
easily shown to have a higher free energy [3l| . and will 
not be considered here. Here and in the following, we 
employ the continuum limit, which also corresponds to 
the infinite bandwidth limit for the electronic dispersion. 
Particle-hole symmetry then implies Ga = 0, whereas 
the fact that the d-wave part of the gap averages to zero 
on the Fermi surface gives G2 = 0. Performing the k 
summation in we find the Green's functions: 



Go(iw„) 
Gi(iw„) 



Al+Lu^ + A? 



-.F{x), (14) 



A, 2No 



A2 + + A2 



-.F{x) (15) 



where Nq is the fermionic density of states, F{x) is the 
complete elliptic integral of the first kind [s^], and we 

defined x = {A\I{A\ + + A^)]^/^ Using the self- 
energy (|10|l . we find that the renormahzed frequency 
is determined by the relation: 



i;(fc, iu^ri) = nmpT{k, k, iun). 



(9) 



For simplicity, we will assume a local potential, v^k' = v, 
so that the solution is easily found: 

A., . , Got" - Girl - Gar^ + {y-^ - G^y^ 

^ ' '"^"^ " ""'"'P -(Go)2 + (Gi)2 + (G2)2 + (G3 ~ 1.-1)2 

(10) 

where we have introduced the local Green's function: 



Gi{iuJn) 



k 



gi{k,iuj„). 



(11) 



Equations |2Jl , ® , COI) > CD constitute the starting point 
of our analysis and must be solved in a fully self- 
consistent manner. In the following two paragraphs, we 
will specialize to the s + idx'^_y2 and d^y + idx2_y2 cases 
respectively, which will allow to carry out further the 
analysis. 



B. Disorder in a s + id^2_y2 superconductor 

We consider here an attractive pairing potential in 
both s and dx2-.y2-'wave channels 



Vkk' 



1.2 _ l2 jL/2 _ i.r2 

f^x + "-J/ '^x + '^y 
Vs + VdCOs(20) COS(20') 



(12) 



F{x 



l + (2A^o«)2(C32 



A2 



F{x) 



(16) 



The renormahzed s-wave gap As = As + Si(ict;„) is given 
by an equation similar to Hlt)|) , which we can write in the 
following way: 



(17) 



This important relation will be central to the later phys- 
ical discussion. Finally, we note that G2 = implies 
S2 = 0, so that the d-wave gap does not acquire a self- 
energy correction in Eq. 

The momentum summation in the gap equation (j^J 
can be performed similarly, and leads to 



1 

7 



E 



2iVo V, As 



-.F{x) 



(18) 



A2+52 



A2 



in the s-wave channel, while the dx2_y2 gap is determined 
by: 



1 ' 

Ad = -- ^ 2iVo Vd Aa 



-w2 + A 



2 r 



A? 



E{x) 



ii2+A2 



-Fix) 



(19) 
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where E{x) is complete elliptic integral of the second kind 

M 

The primed summation indicates that a frequency cut- 
off ojc needs to be introduced (see Appendix . We 
emphasize that this procedure is crucial in order to pre- 
serve Anderson's theorem in a superconductor with non- 
retarded interaction such as (|12|l . Indeed, the use of a 
finite energy bandwidth together with an unbounded fre- 
quency sum would violate this basic result 0, • Al- 
ternatively, one could employ a retarded pairing interac- 
tion from the outset, without restricting the frequency 
summation, however, this procedure would significantly 
complicate the analysis without modifying the results. 



Disorder in a d^y + id^ 



2 superconductor 



We now consider an attractive pairing potential in both 



dxy and dx2_y2--weive channels: 

Vkk' — Vxy 1 r^if y ~^ 

^l%xf^y ^f^x'^y 



Vd 



h.2 I 1,2 i,/2 I 1.12 



^(20) 



= Va:y sin(20) sin(2(/)') + Vd cos(20) cos(2(?!)'). 

This leads to a d^y + id^2_y2 superconducting order: 

Afc = A^y sin(2(/)) + iAd cos(20). (21) 

Arguments discussed previously give Gi ^ G2 = G3 = 0, 
so that both A^y and are not renormalized. The only 
non-zero local propagator is: 



Go{iUn) 



iu 2Nq 



2 

M 



fFiy) 



(22) 



'2M-1/2 



where y = [(A|, - A^ ^/(A^^ + ^^)] with the def- 
initions Am = M(ix{A^y,Ad) and A„ = Miii{A^y, Ad). 
There is therefore a single equation to solve for the renor- 
malized frequency: 



Fjy) 



1 + {2Novyi:j 



2-2 F(y)^ 



(23) 



D. Results for the s + id^2^y2 case 

We start by reviewing the occurence of the quantum 
phase transition in a clean superconductor with two com- 
peting pairing channels. Then we study the fate of the 
critical point in the presence of disorder. Finally we 
present the phase diagram obtained from the numerical 
solution. 



1. Clean limit 

It is useful to recall the nature of the ground state in 
a superconductor with competing s and d-wave instabil- 
ities. Let us first consider the possible occurence of an 
s-wave instability in a clean well-formed (ij,2_j,2-wave su- 
perconductor. Basically, the BCS singularity is cut off 
by the presence of the dx2_y2-wave gap, as shown by the 
fact that the propagator behaves as {A'^+u;'^)~^/^ in the 
gap equation H18|l . Putting As to zero in this equation, 
one sees that a finite critical interaction Vs is necessary 
to trigger the s-wave instability, in contrary to the situ- 
ation in a normal metal. From the same argument, the 
s-wave state is stable at small Vd up to a critical value of 
Vd where the additional d-wave order parameter devel- 
ops. It is straightforward to establish that there exists a 
finite range domain in which both orders coexist at zero 
temperature (see (31] and Appendix : 



\Vd\ 



2 + No\Vd\/2 



< \Vs\ < (26) 



For instance, the quantum critical point associated to the 
transition dx2^y2 — > s + idx2-y2 occurs when couplings 
obey \Vd\/{2+No\Vd\/2) = \Vs\ and can be described by a 
field theory involving the coupling of nodal quasiparticles 
to Ising fluctuations associated to the additional s-wave 
order ^| • - Clearly, this description lies beyond the 
scope of the present mean-field analysis. 



2. s-wave instability in a dirty d^ 



-wave background 



Setting Vm = Max(Kj,,yrf) and V,n = Min{V^y,Vd), 
we finally express the two gap equations: 



A„ 



lY.2N,Vrr.A^^ 



M 



Al 



Am = II 2NoVmAm'^ 



A2 



A^ 



F{y) - E{y) 
(24) 

E{y) 

(25) 



We now perform the analytical and numerical analysis 
of the previous BCS-|-T-matrix equations, considering in 
turn the s -|- idx2_y2 and d^y + id,j.2_y2 cases. 



Let us start with an important remark on the effect of 
disorder in s-wave superconductor {Vd = 0). The local 
propagator ifT^ reads 



Go{iLUn) = - 



iuj t:Nq 



iuj ttNq 



32 , A2 y^^TAf ' 



(27) 



where we have used the relation (|17() to get the r.h.s. 
expression. All renormalization factors in Gq drop out, 
and therefore the gap equation p8|l of the pure s-wave 
superconductor is unchanged by the presence of disorder, 
in accordance with Anderson's statement. For the pure 
d-wave case, no such cancellation occurs, and the gap 
is rapidly suppressed when a few percent of impurities 
are added. We note that, in reality, Anderson's theorem 
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holds to order nimp [3^ , so that the strict insensitivity of 
As to disorder here is actually an artefact of the T-matrix 
approximation. 

We now consider a c?-wave superconductor (V, — 0) 
containing a small amount of disorder (such that the gap 
Ad is not completely driven to zero). From (|16|l . assum- 
ing unitary scattering (w = oo) but small impurity concen- 
tration (riimp ^ 1), we find the renormalized frequency 
uj = ^ nimpAd/(2A^o) at = 0^ (up to logarithmic cor- 
rections coming from the elliptic function) . This leads to 
the density of states at the Fermi level: 



po = --ImGo(iO+) 

TT 



2A^o 



2A (28) 

i.e. a finite density of states is generated at low energy 
due to the pair breaking effect of the non-magnetic impu- 
rities [s^lSq . This result is also valid in the Born limit of 
small scattering potential u, although po is exponentially 
small in riimp in this case j22j . Experimentally, scatterers 
in high- Tc materials appear to be closer to the unitary 
limit. 

Now, we focus on the s-wave attraction and study 
the related BCS equation (|18|) in the presence of this 
disorder- induced density of states po- By similarity to 
the BCS instability in a normal metal, we would naively 
say that any infinitesimal would lead to a finite s- 
wave gap A^. As we will see in the later case of the 
dxy + id.j;2_y2 superconductor, this argument is true pro- 
vided Anderson's theorem is fulfilled, which is the case for 



superconductors only. Indeed the BCS equa- 



tion H18|) can be re-expressed using again relation (|17|l : 



2jVo Asu/lu 

/3 ^ ^AjT^^TApv^ 



2No Vs Asw 



■Fix) 



:F{x). (29) 



/3 ^ yZJ5(A2T^5)TAp 

Because uj is finite at zero frequency in a dirty d-wave 
superconductor, a l/w singularity is generated at van- 
ishing As and the weak-coupling BCS instability to the 
s + id^2_y2 state occurs, in accordance with our intuitive 
argument. Let us summarize this interesting result: due 
to Anderson's theorem, an arbitrarily small s-wave com- 
ponent of the pairing potential will always generate an 
s-wave contribution to the order parameter in a disor- 
dered d-wave superconductor. Of course the critical tem- 
perature for the onset of the s-wave gap is exponentially 
small in the disorder-induced density of states (up to mul- 
tiplicative factors): 



s— wave 



■ exp 



/2AdiVo 




(30) 



-ujc exp 



TT^Arf 



2No\Vs\n 



imp 



In conclusion, this disorder-induced s-wave instability 
could be observed in practice by increasing the impu- 
rity concentration riimp. Note that, while the exponential 




0.08 



FIG. 1: (color online). Superconducting gaps Ad (continous 
lines) and As (broken lines) as a function of nimp for various 
values, from left to right, of the s-wave pairing Vs = —0.45 
(mixed region), Vs = —0.40 (s ^ s + id critical point), Vs = 
—0.3 (clean d phase with s component induced by disorder) 
and Vs = —0.1 (pure d phase). Here Vd = —1.0 is fixed and 
the inverse temperature is taken as /3 = 1000. 



term in 1)301) raises from zero by introducing disorder in 
the system, the global prefactor is suppressed by the y/Ad 
term, and T^~""'^'^^ will be maximal for an intermediate 
value of n;,nn. 



3. Phase diagram 

We present now the numerical solution of the self- 
consistent BCS-I- T-matrix equations. In all following 
calculations, = 1 is chosen as basic unit, and the 
frequency cutoff in the Matsubara sums is ojc — 0-1 
(such a low value of ujc allows to obtain realistic gap 
estimates 0). The unitary limit [v = oo) is also con- 
sidered here. Importantly, the physics discussed below 
will be qualitatively insensitive to specific values of these 
parameters. Considering the mean-field nature of our 
model we will not try to make quantitative contact to 
experimental data. 

Fig-Dillustrates our main result by plotting the evolu- 
tion of Ad and As as a function of impurity concentration 
J^imp- We notice that A^ is generically suppressed by the 
addition of disorder, whereas As is generated even if not 
initially present in the clean system (assuming 7^ 0). 
Because these curves are computed at finite tempera- 
ture, a critical interaction Vs is necessary to enter the 
s -|- id^2_y2 state, as shown in the related phase diagram. 
Fig. 13 However, when temperature is lowered, the pure 
d3;2„j,2-wave phase is gradually suppressed due to the un- 
avoidable s-wave instability discussed previously, leading 
to the interesting zero-temperature phase diagram shown 
in Fig.El 
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FIG. 2: (color online). Phase diagram (|Vs|,ninip) of the 
dirty s + id^2_y2 superconductor for (3 = 1000. Thin dashed 
lines refer to the gaps plotted in Fig. The inset illustrates 
the evolution of the phase boudaries with temperature, here 
shown for /3 — 100 (broken lines) and (3 — 1000 (continuous 
lines) . 



FIG. 4: Superconducting gaps (continous line) and A^y 
(broken line) as a function of Wimp for the dxy-WtWG pairing 
Vxy = —0.47 (mixed phase in the clean limit), with Vd = —0.5 
and (3 = 100. 




dx^-y^ phase obey Anderson's theorem. This will have 
important consequences for the following analysis. 



2. dxy-wave instability in a dirty d^2_y2-wave background 

We now consider the effect of a small pairing ampli- 
tude in the dxy channel, supposing that the dx2-y2 order 
is present despite the finite amount of disorder. Although 
a finite density of states po is present, we easily see that 
the naive argument predicting a weak-coupling BCS in- 
stability for the /S.xy gap does not hold. Indeed the gap 
equation H24|) in the dxy channel reads: 



FIG. 3; (color online). Zero-temperature phase diagram 
(|Vs|,nimp) in the dirty s + id^2_y2 case. The thick line stands 
for the limited pure d^2_y2 region. Compare with the finite 
temperature plot displayed in the inset of Fig. |5| 



E. Results for the dxy +id^2_y2 case 

1. Clean limit 

In the absence of disorder, dxy and dx2-y2 order pa- 
rameters compete as in the previous case, leading to a 
pure dx2-y2--weLYe state at small Vxy, a pure dxy-'wetve 
state at large Vxy and a mixed region of dxy + idx^-y^ 
symmetry in between (the dxy + dx-2-y2 state being un- 
favorable energetically). This region is found to be given 
by the condition (see Appendix IXjl : 



(31) 



Of course, neither the pure dxy phase nor the pure 



^xy^-^Yl Vxy Ax 



A2 



Fiy) 



E{y) ■ 

(32) 

The integral is completely regular at low frequency, and 
Axy can be put to zero, leading to a critical interaction 
Vxy to generate the c?a;j,-wave gap, similarly to the clean 
case. The physics is therefore very different from the s + 



idx 



„2 situation, and is due to the absence of Anderson's 



theorem in dxy-^ave superconductors. 



Phase diagram 



The previous discussion implies that the results for the 
y2 case are markedly different from the one 

nci c _1_ ^ n ^ 



dxy + id.. 



obtained in the s -\- idx^-y^ case. Indeed, the numerical 
solution of the mean-field equations show that, if the dxy- 
wave gap is not formed in the clean limit, it cannot be 
generated by adding impurities. By the same token, in 
the mixed region, both order parameters and Axy are 
gradually suppressed. This is shown in Fig. ^ 
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FIG. 6: Gap structure of a bond-ordered superconductor 
(dubbed spin Peierls in the following). In general, gaps take 
complex values, so that time reversal symmetry is also broken. 



FIG. 5: Phase diagram (|'l4y|,7iimp) of the dirty dxy + idx2_y2 
superconductor, with Vd = —0.5 and (3 — 100. The thin 
dashed line refers to the gaps drawn in Fig. |1| 

The phase diagram illustrates this behavior, see Fig.|Sl 
We stress that the zero-temperature phase diagram is 
quite similar to this picture, i.e., the pure dx2_y2-weLve 
region remains extended, in contrast to the s + id^2_j,2 
case. 



III. C-BREAKING IN DISORDERED 
SUPERCONDUCTORS 

In contrast to conventional superconductors, cuprates 
derive from their Mott insulating parent compounds, and 
the strong-correlation physics of doped Mott insulators 
is responsible for many of the intriguing phenomena ob- 
served in these high- Tc materials. Theoretical studies 
of microscopic models, such as the t — J and Hubbard 
models on a two-dimensional square lattice, have shown 
that a variety of ordering tendencies compete, and that 
the actual ground state may depend sensitively on micro- 
scopic details. Experimentally, some compounds shown 
evidence for lattice symmetry breaking in form of charge 
and/or spin order, possibly coexisting with superconduc- 
tivity at low temperatures. 

In this section, we shall consider the competition be- 
tween superconductivity and additional order associated 
with broken translation symmetry under the influence 
of disorder. Many candidates for the additional order 
could be chosen [l^ , but we will concentrate here on the 
simplest case of additional bond (or spin-Peierls) order, 
where bond strengths are modulated as shown in Fig. 
Such a state obeys a two-site unit cell, but up to a rota- 
tion all sites are equivalent. This justifies the application 
of the T-matrix approximation, based on the exact so- 
lution of a single defect, where the use of a single self 
energy is sufficient. In contrast, for more complicated 
forms of order, such as stripe states with varying site 
charge density, the T-matrix approach would clearly be 
inappropriate, and this is one of the reasons to restrict 



our attention to bond order. 



A. Formalism 

We start by developing the necessary formalism in the 
absence of disorder, then compute the T-matrix in such 
a non-translation-invariant superconducting background, 
which will finally lead us to study how disorder affects the 
chosen ordered state. 



1. Clean case 

Contrarily to Sec. ^ we cannot directly simplify the 
analysis in the continuum limit as the lattice symmetry 
is now explicitely broken. The simplest discrete model 
that exhibits the physics we are interested in is given by 
the t — J Hamiltonian: 

RR'(T 

(33) 

where the operators prohibit from double occu- 

pancies, the rest of the notation being standard. Re- 
stricting both hopping t and exchange interaction J to 
nearest-neighbor sites, and introducing a slave boson 
6^ to implement the non-double occupancy restriction, 
with the representation = /r^-^r ^-nd the constraint 
Y.a /L/r<t + ^R^R = 1, we arrive at 

H = -t ML'/L/r'.-mE/L/r. (34) 

(RR')(T Ro- 
(RR')cts 

where (RR') denotes nearest neighbor pairs of lattice 
sites, and fi is the chemical potential enforcing the con- 
straint. We have organized the interaction term in such a 
way that the mean-field treatment in the particle-particle 
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sector (controlled by a large- limit in the Sp(7V) formu- 
lation 13,^3) is quite straightforward: 

Hmf = -toff E fLfn'.-^^T.fLfIL. (35) 

(RR')o- Ro- 
(RR')o- 

with 

Arr, = -JE(-1)'(/r'-s/rs) (36) 

s 

This mean-field is equivalent to the pairing part within 
the BCS formalism (which can be argued to be 
enough even when pairing is induced by the presence of 
strong correlations). The slave bosons 6r in (|34() con- 
dense at low temperatures, and at the mean-field level 
the 6r operators are simply replaced by their conden- 
sation amplitude 6^ = (&r) = v^, where S is the hole 



doping, J2Ra f-Rcrfna = 1-^- Thus, the strong Coulomb 
repulsion in the model is accounted for by a renor- 
malized hopping tcs — tS in l|35|l . 

If we assume the symmetry breaking pattern shown in 
Fig. El we have two sites per unit cell. On these two 
sublattices, we can define the fermions = c\j^ for 

R = -I- jey and /^^ ee d\.^ for R = (2/ + l)e^ + 
jBy, where I and j are integers. Generalizing the Nambu 
formalism, we introduce a four-component spinor: 

with Kx = 2k,j: the supermomentum associated to / 
and ky the usual momentum associated to j. After 
some manipulations (see Appendix^, we find the mean- 
field Hamiltonian in matrix form Hmf = X^if k ^^^fc^i 
where: 



hh = 



— 2toff cos ky — ji 
2A3 cos ky 



■26 



2 A3 cos ky -tcs - tcse 
2tcs cos ky + n Ai+ AJe" 



Ai 

teS 



Ai 



^26 



-2tcscosky - fi 

2A3 COS ky 2teS COS ky + fi 



2 A3 COS ky 



(38) 



The gap equation 136|) can similarly expressed (see Ap- 
pendix in terms of Green's functions of the spinor: 

= iujn - /ife - E (39) 

where E contains effects due to the disorder distribution, 
that we will evaluate now. In principle, E has also con- 
tributions due to fluctuation effects beyond mean-field, 
which we will not take into account here. 



2. Calculation of the T-matrix 



3. Self- consistency cycle 

In the presence of a finite but small concentration 
n-imp of impurities, we can approximate the self-energy by 
E ~ ninipT. Self-consistency must be reached at two lev- 
els: first because the T-matrix is expressed through the 
propagator by the result H4U|I and Dyson's equation H39|l . 
second because the gaps Ai, A2, A3 should obey the BCS 
equations (|B3IB5ll . This whole procedure is performed 
numerically, with a discrete momentum mesh and at 
small finite temperature. 



The calculation of the T-matrix corresponds to solving 
the problem of a single localized impurity at R = in the 
superconducting background state. It is obtained from 
the solution of the problem Hmf + "^^a foafoa' 
after some manipulations, we find (see Appendix IB|) : 



1 



GiiG 



22 



G12G 



12'-J21 



~G22 


G21 








G\2 


Gil 














— G44 


G43 








G34 


G33 



(40) 

where Gm„(ia;„) = Y^K^ky 9"inii^n, Kx, ky) is the local 
propagator. For simplicity, we have assumed the unitary 
limit V = +00. 



B. Numerical results 



Clean case 



In the absence of disorder, the Hamiltonian H34|) shows 
several phases as a function of the pairing attraction J 
and doping J H [H. The C-broken bond-order phase, 
which coexists with superconductivity, occurs typically 
at low doping, and is characterized by the Ai, A2, A3 all 
being different. Because the effective hopping t^ff = tS 
can become much smaller than J in this regime, this 
corresponds to a strong-coupling superconductor. The 
bond-ordered state also breaks T, as the three bond vari- 
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FIG. 7: Superconducting phase diagram of the clean t — J 
model for /3 = 100 (in units of J = 1) , with a two-site unit cell. 
The transition from s* to d^2_y2 is discontinuous, whereas all 
other phase boundaries are second order. These results are 
consistent with 

ables cannot be made all real by a global gauge trans- 
formation. At intermediate doping, the ground state is 
uniform, but continues to break time-reversal symmetry, 
as the s* + id^2_y2 phase happens to be more stable (if 
J is large enough only, otherwise the pure (i^2_j^2-wave 
state with Ai = A2 = —A3 is found). Finally, an ex- 
tended s-wave phase, denoted s* and characterized by 
Ai = A2 = A3, is encountered at larger doping. We have 
checked numerically these known results for the clean 
case, as shown in Fig. [7| 

2. Dirty case 

We investigate now how these phases evolve with the 
addition of non-magnetic impurities. For this model cal- 
culation, we stress that Anderson's theorem does not 
hold, as a finite energy bandwidth is used (see the dis- 
cussion at the end of Sec. Ill B|) . For this reason, and in 
contrast to the results obtained in the first part of the pa- 
per, we do not see transitions from d^2_y2 to s* -\-id^2_y2 
by increasing the disorder, and rather all superconduct- 
ing phases are destroyed in a similar fashion. 

Let us discuss the effect of adding impurities to the 
bond-ordered phase. The numerical solution of the prob- 
lem demonstrates that the translation symmetry is re- 
stored during this process, i.e. the C-broken phase is 
generically unstable to the addition of disorder, as dis- 
played in Fig.|Hl This result can be expected on the basis 
that impurities disrupt the perfect arrangement of super- 
conducting dimers, suppressing their long-range order. 
(In principle, it is possible that dimers locally arrange 
around i mpu rities, resulting in enhanced local dimer cor- 
relations 1,221, but this is not related to bond long-range 
order, and is certainly beyond the scope of our mean- 
field approach.) Discussing in more detail the evolu- 




0.1 0.2 0.3 



FIG. 8: (color online). Superconducting gaps as a func- 
tion of impurity concentration, for t/J = 1, /3 = 100 and 
S — 0.375. Single (double) primed quantities denote real 
(imaginary) parts, respectively. Here we chose a gauge where 
A'l = 0. The initial state at nimp = is the bond-ordered 
superconductor described in the text. 
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^imp 



FIG. 9: (color online). Superconducting gaps as a function 
of impurity concentration, same conditions as in Fig. m with 
5 = 0.45 here. The initial state at nimp = is s* -f id^2_y2. 



tion of the bond parameters shown in this figure, we see 
that the strongly dimerized superconductor persists for 
< nimp < 0.07, while being progressively suppressed. 
We then have a narrow s* -\- id^2_y2 phase with restored 
lattice symmetry for 0.07 < riimp < 0.10 (in which 
A'/ = A^' = and A^ = A'^ = lA'^ + iA'^\). Finally, a 
d3;2_j,2-wave state (gaps are real, with Ai = A2 = —A3) 
is found at larger impurity concentration, ultimately un- 
stable towards the s* phase (with Ai = A2 = A3) for 
njinp > 0.25 in a discontinuous manner. 

Similarly, starting at larger doping from the clean 
s* + id^2_y2 phase. Fig. |^ shows an initial suppression 
with riiiiip of the T-breaking towards the dx2-y2-wsLve 
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FIG. 10: Superconducting phase diagram of the dirty t ~ 
J model as a function of impurity concentration riimp and 
doping 5, for t/J = 1 and (3 = 100. The two dashed lines 
refer to the cuts in Figs. |H] and |5] respectively. 



phase, which finally displays a first-order transition into 
the extended s*-wave phase. In all these calculations, 
we note that the unrealistic value for the critical concen- 
tration, where all superconductivity ultimately vanishes, 
is an artifact of the mean-field approximation: to stabi- 
lize the bond-ordered phase, a very small t^g is required, 
placing the superconductor effectively into the strong- 
coupling limit which is rather robust against impurities. 
Nevertheless, we expect our results to correctly capture 
the qualitative aspects of impurity doping. 

The described sequence of phases as a function of im- 
purity concentration can roughly be guessed from the 
knowledge of the phase diagram in the clean case. We 
summarize our results for the dirty C-breaking supercon- 
ductor in Fig. ^1 We have not mapped out the full set 
of phase boundaries for all possible parameter combina- 
tions, thus for different values of t/J, the various phases 
might organize themselves in a different way. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper, we have investigated the influence of 
static disorder on superconductors with competing in- 
stabilities. Regarding the possibility of T-breaking states 
with a mixed gap of s -I- id^2_y2 and dxy + idx2_y2 struc- 
ture, we found important qualitative differences between 
those two cases in the process of adding non magnetic im- 
purities. In studying next a superconductor with bond 
long-range order (doped spin-Peierls state), we found 
that the lattice symmetry is very quickly recovered by 
addition of disorder, and we would expect this effect to 
be generic. 



A. Application to high-temperature 
superconductors 



The finite-temperature phase diagrams displayed in 
Fig. [3 for the s + idx2_y2 superconductor and Fig.|5lfor 
the dxy + idx2-y2 superconductor show clear qualitative 
differences to the effect of impurity addition. Indeed, the 
minor s-component is enhanced with respect to the dom- 
inant dx2 -y2-wave background in the first case, whereas 
the small d^y part of the gap is very rapidly suppressed in 
the second case. This implies that typical tunelling spec- 
tra would evolve in a very different manner with impurity 
concentration, depending of the underlying structure of 
the superconducting gap: in the s + ic?^2_j,2 case, we ex- 
pect that the zero-bias conductance peak splitting should 
increase, while the opposite behavior could be seen in the 
dxy+idx2-y2 case. It would thus be very interesting to see 
such an experimental study in overdoped YBa2Cu307_5, 
which could lead to a definite statement on the nature of 
the mixed superconducting gap in this compound. We 
recognize, however, that the practical situation might 
be complicated by the fact that magnetic moments can 
be induced by non-magnetic scatterers such as Zn [38j |. 
which then would strongly affect the s-wave component 
as well. 

A second point we would like to emphasize is that 
the typical coexistence windows given in equations l|26() 
and H31(l are extremely narrow in the weak-coupling 
regime NoVd <C 1. Because of their smaller critical tem- 
perature, we believe that electron-doped superconductors 
as Pr2_a;Cea;Cu04 fall into this condition, explaining the 
absence of clear observation of a mixed phase in tunnel- 
ing experiments [15| . 

Let us turn to states with additional breaking of trans- 
lation symmetry, i.e., charge order. Our calculation 
nicely shows that the global symmetry is restored very 
quickly upon adding impurities. Clearly, physics beyond 
mean-field can lead to the nucleation of local symmetry- 
breaking patterns around impurities due to pinning ef- 
fects. Then, in the presence of a sufficient amount of 
impurity disorder, signatures of charge ordering will only 
be visible in local probes, such as STM, but no sharp 
superlattice Bragg peaks will be seen, and no thermo- 
dynamic phase transition will occur. This scenario defi- 
nitely has similarities with what is seen in certain cuprate 
compounds [2lll2ll2^. 

Competing orders are ubiquitous to many exotic su- 
perconductors. Thus we think the theory developed in 
this paper would likely find applications in some other 
context. Candidate examples are Naa;Co02 [13 , a pos- 
sible realization of a doped Mott insulator on the trian- 
gular lattice, and organic superconductors of the BEDT- 
TTF family [43| , where both spin and charge fluctuations 
seems to play a crucial role for the pairing mechanism. 



11 



Acknowledgments 

We thank P. Hirschfeld, X. Wan, and P. Wolfle for 
helpful discussions. This research was supported by the 
Deutsche Forschungsgemeinschaft through the Center for 
Functional Nanostructures Karlsruhe. 



APPENDIX A: DERIVATION OF THE GAP 
EQUATION FOR T-BREAKING 
SUPERCONDUCTORS 



Performing the oj-integral in the large ujc limit, we find: 

V-^~2V^^ = -2iVo / d<p cos(20) log (a^ + cos^ 0) 
Jo 

= ^ (^-l-2a- + 2aVa^ + l) (A3) 

where a = A^/Arf. The instability criterion for the s- 
wave phase (resp. dx^-y^-^SLve) is given by a = (resp. 
a = oo), leading to the mixed region H26(l . 



1. s + id^2_„2 superconductor 



2. dxy + id^2_y2 superconductor 



In order to set up clearly the notation, we present how 
the gap equations H18|) - (|19|l can be obtained from 
We first transform the k sums into an integral over an- 
gle and energy in the infinite bandwidth limit, allowing 
Anderson's theorem to be preserved |3J, |35| : 



A. 



A. - - 



(3 ^ 

n 

-e' 



27r 
2^ 



+ 00 

de- 



-oo w2 + A2 + A2 cos2(2<?!)) + e2 
de- 



NpVdAd cos^(20) 
A2 + A2cos2(2</))+e2 



Because of the self-consistency, uj and Ag are unknown 
functions of the frequency a;„, and the Matsubara sum 
cannot be performed analytically. However, the integral 
over e is straightforward to do: 



A,, = - 



-E 



tt/2 



2NoVsAs 



A. 



/3 



E 



7r/2 



j2 + A2 + A^ cos2 , 
2NoVdAdCos^(t) 



A2 + Aj cos2 , 



(Al) 



(A2) 



which leads indeed to expressions l(TI^ - ((T^ . 

The 1/uj large frequency behavior in the Matsubara 
sum is due to taking the infinite bandwidth limit, and 
a frequency cutoff ujc must be introduced. In practice 
we have not employed a hard cutoff, as this leads to a 
first-order jump of the order parameter at the critical 
temperature, but rather a smooth cutoff function, 
J2nfi^n), with: 



1 if LJ < UJc 

g-2(t^-t^e)/t^, if W > 



The derivation of the gap equation is very similar to 
the previous case, and we only quote the intermediate 
result: 



Axy — 



A, 



7r/2 



2N(,VxyAxy sin^ 



w2 -I- A2 , sin^ + A^ cos^ , 



(A4) 



E 



d(f> 





2Aro^dA<jcos2 



a;2 4- Aly sin^ (p + cos^ ( 



(A5) 



which gives expressions (j24|l - (j25|l after some manipula- 
tions. 

We also find the location of the critical points in the 
clean limit. Subtracting l)A5|l from (|A4p at zero tem- 
perature and performing the w-integral at large uJc leads 
to: 

V-y' - V^-'^~2NoJ^ d</)cos(2<^) log (^Y3^ + co^' ^) 

(A6) 



2 a+1 

where a = A^y/Ad now. The instability criterion for the 
dxy-'wsLve phase (resp. dj:2_y2-wave) is given by a = 
(resp. a = oo), leading to the mixed region H31|) . 



APPENDIX B: DIRTY SUPERCONDUCTOR 
WITH BOND ORDER 

1. Gap equation 



Finally, let us determine the location of the critical 
points in the clean limit 'si']. Subtracting twice (|A2|I 
from IjAljl at zero temperature gives: 



V-^-2Vi^ = 2No 



r/2 



We give here some details on the derivation of the gap 
equation for the bond-ordered superconductor. Due to 
the broken lattice symmetry, the resulting 2x1 unit cell 
leads us to introduce two fermion species on each sub- 
lattice, c\j^ and d\j^, and we can directly express the 
27r ^Ci;2 -\- A'^ + A'^ cos^ (j) Hamiltonian (|^ in real space in those variables. It is 



duj 



cos(2(/)) 
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then convenient to go to Fourier space, and we set 

4. = EE4..„.e^^-^^+'=«^ (Bl) 



'■Ija- 



EE4 



KxkyCr 



(B2) 



where = 27r/(L/2), . . . , 27r and ky — 27r/L, . . . , 27r 
take respectively L/2 and L values (L x L is the total 
number of sites). Inserting these expression in the mean- 
field Hamiltonian gives readily the final matrix expres- 
sion 

We can also write the gap equation H36|l in terms of 
these degrees of freedom: 

^1 ^ ~'p E "^^^^y i'^-K^-ky-a'^K^kya) 



A, 



A, 



J2 Y.(~ir^"'^(' 



-K^—ky — a Kxk 



iy (7 } 



, ky a 



Introducing the propagators gmn = (^n^'m): we can 
write the gap equations as: 



Ai = E [514 + 532] 



(B3) 



LJn Kxky 



purity at site R = 0, and start with the exact Dyson's 
equation [iuj - Hmf - vJ2afoafoa]9 = 1' which we 
want to sandwich between two Bloch states \k,n^, where 
k = {Kxj ky) and n = 1, . . . , 4 is the Nambu index. In- 
troducing the propagator (7m„(fc, A:') = (fc, TO|g|fc', n), we 
find in matrix (Nambu) notation: 



[iuj-hk\g{k,k') = 5kk- +vY,Dg{k\k') (B6) 



D 



(B7) 



1 0- 
0-100 

0. 

Performing the matrix algebra, we can solve the previous 
equation, and we find the local propagator: 



G{k,k') 
f{iuj) 



Jkk' 



f{iuj)- 



lu) — hk' 
1 



1 

i5et 



k 

-1 _ 

G12 





iu - hk 



-D 



(B8) 
(B9) 



G 



22 



G 



21 





Gil 







A, 



A, 



J 



J 



E [e'^"5i4 + e '-^==532] (B4) where we have defined: 



^K^ky 



2 COS ky gi2 



(B5) 



^K^ky 



2. T-matrix 



Det = GiiG22-Gi2G2i-i'"^(Gii + G22)- 



We prove here the formula (|40|l for the T-matrix of a 
bond-ordered superconductor. We consider a single im- 



(BIO) 



Taking into account the contribution from site R = e^;, 
and in the unitary limit w = 00, we find the final expres- 
sion igOI). 
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